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We predict that large moments J, placed into a crystal 
field with the cubic point symmetry group, differ by their 
spectrum and magnetic properties. E. g., properties of the 
odd-integer moments are different from those of the even- 
integer. The effect is due to Berry's phases gained by the 
moment, when it tunnels between minima of the external field. 
Two cases of the group O are classified, namely, 6- and 8- 
fold coordinations. The spectrum and degeneration of energy 
levels depend on a remainder {J/n}, where the divisor n — 4 
and 3 for 6-fold and 8-fold coordination respectively. Large 
moments in the cubic environment can be realized by diluted 
alloys Ri-xR^Sb, where R=Lu, La, and R'=Tb, Dy, Ho, Er. 

PACS numbers: 03.65.Bz, 75.30.Cr, 75.10.Dg 



It has been shown on several examples (dependence of 
the Haldane gap in anti-ferromagnetic spin chain Q and 
in 2D isotropic Heisenberg antiferromagnets (2j , suppres- 
sion of tunneling in magnetic particles |^] , magnetization 
of ions in external fields [Q ) that systems with large inte- 
ger moments may behave differently from ones with large 
half-integer moments although the values of the moments 
are close to each other. Haldane pointed out in g] that 
there may be a difference between odd and even spins as 
well. All of them are due to Berry's connection ||. In 
this paper we demonstrate systems where further differ- 
entiation on the moment values occurs. 

The systems under study are large moments in strongly 
anisotropic external fields, e.g, crystal electric field with 
deep wells at certain crystallographic directions. We con- 
sider systems of the octahedron symmetry with 6- and 
8-fold coordination. They show a subtle partition on the 
moment values. 

In classical description, the angular moment is a vec- 
tor of a fixed length J determined by the direction of a 
unit vector n. Let us assume that an ion with a large mo- 
ment is placed into a crystalline electric field (CEF). Due 
to symmetry the moment may be localized in 6 symmet- 
rical positions (the easy positions of CEF). We choose 
following spherical coordinates (6, 4>) for the positions: 
l-(0, 0), 2-(tt, 0), 3-(7r/2, 0), 4-(7r/2,7r), 5-(vr/2, tt/2), 
and 6-(7r/2, 37r/2). Such a configuration will be called 
octahedral or 6-fold. The ground state of the system is 
6-fold degenerate if tunneling between the potential wells 
is neglected. If one allows the moment to tunnel to the 
nearest wells along the geodesies on the unit sphere (see 
Fig. [l]) , the degeneracy will be lifted at least partially. 




FIG. 1. Paths of the spin on the unit sphere between the 
easy positions of the field. The case of the 6-fold coordination. 

A formal description of such a system is given by a 
Hamiltonian (H) with the following matrix elements: 



hi 
hi 



0, i = 
0, It 



1, 2, ...6, 



3, 7, 11, 



(1) 



| w| , for other 1 < i,j < 6 



where we adopted the enumeration shown on Fig. [l| 

The moment's wave-function, upon tunneling from one 
minimum (i) to another (J), may gain a geometric phase 
shift <p(Cij) which depends on the tunneling path Cy 
on the unit sphere (the n-space): 4>(Cij) = J c A(n)dn. 
The vector field A is called Berry's connection. Specially 
for the problem of semi-classical moments it was calcu- 
lated in . The appropriate geometric phase was first 
discovered by Berry (|] . 

The phase shifts 4>{Cij) as well as the connection A 
are not gauge invariant: A(n) — * A(n) — V/(n), where 
/(n) is an arbitrary continuous function, the spins' eigen- 
functions transform as: ^(n) — ► e*-^ n ^(n). In our case 
(well localized states) it reduces to: Vfc — * e^ fc Vfcj where 
fk = /(n/-), rife (fc = 1, ... ,6) are the positions of the 
potential wells. However, the phase 7(C) accumulating 
at circulation along a closed circuit C is invariant: 



1(C) 



j> A(n) 



dn = ±m(C), 



(2) 



where f2(C) is the solid angle subtended by a closed path 
C on the unit sphere at its origin. The phase shifts <fi(Cij) 
change to <p(Cij) + fj — fi by the gauge transformation. 
Therefore, 4>(Cij) can be any set of numbers satisfying 
the invariant relationships (|^). The Hamiltonian also is 
not invariant, it transforms as: "H — > WHU\ where U is 
a diagonal matrix with matrix elements Ukk = exp(i//-). 

Tunneling trajectories divide the sphere into 8 pla- 
quettes. Each of them represents a minimal non-trivial 
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closed loop and subtends the solid angle of 7r/2 (e. g. 
closed path 1-3-5-1 on Fig. 0). The 8 plaquettes give o 
equations for the phase shifts: 

4>ij + 4>jk + 4>ki = J{-z mod47r), (3) 

where tf>ij is the phase shift upon tunneling from easy 
position i to easy position j, and set {(i,j,k)} enumer- 
ates the plaquettes. Only 7 out of the 8 equations are 
independent which leads to 5 degrees of freedom over 12 
phases. This number (five) is equal to the number of 
independent parameters in a general gauge transforma- 
tion (six) minus one corresponding to a common phase 
factor. Note that the Hamiltonian ([!]) is not invariant 
under the action of the octahedron group realized as a 
subgroup of the permutation group on the 6 states. What 
leaves the Hamiltonian invariant is a rotation (TV) ac- 
companied by an appropriate gauge transformation (U): 
Ti = U1ZH1Z T U^ , where U depends on TZ, and J. 

Hamiltonian (Q) can be diagonalized leading to the fol- 
lowing 6 eigenvalues: 

E k (J) = (-l) k 2w X (n(J + 2k)), k = 0, . . . , 5, (4) 

9 2x , i\ 

x) = cos — cos cos — h sin — sin — 

' 3 2 V 3 3 2 J 

The ground state energies and their degeneracies fall into 
5 classes which are gathered in Table Bj. The classes dif- 
fer by their ground-state energies, degeneracies of the 
ground-state, and excitation energies. All these will lead 
to experimentally observable dissimilarities between sys- 
tems falling into the different classes for k^,T < w, e. g, 
thermodynamical properties or response to applied mag- 
netic field. 

The number of the classes and the moments' values of 
each class can be easily found without detailed analysis 
of the Hamiltonian for the external fields of the perfect 
polyhedra group symmetries (perfect Hamiltonians) . Let 
us assume that the parameter space (unit sphere) can 
be covered with p plaquettes (the number of the faces 
of a perfect polyhedron). Then, each plaquet subtends 
a solid angle of Air/p at the origin of the sphere and 
Berry's phase for each loop is AirJ/p. This gives a pe- 
riod of J: J p = p/2. The spectra of the systems differ- 
ing by transformation 7(C) — > —7(C) must be identical 
due to time-reversal symmetry. It leads to equivalence of 
J = JpTi + m/2 and J = | J p n — m/2\, n, to = 0, 1, 2, . . .. 
Hence, we can classify equivalent moments in the follow- 
ing manner: 

Tfl 

J m = \ J P n ± — I, n = 0, 1, 2, .. .; m = 0, 1, . . . , J p , 

and the number of the classes is C p = J p + 1 = p/2 + 1. 
The classes can be conveniently labeled by the integer 
to + 1. E. g.. for the 6-fold coordination (octahedron) 
one has: p = 8, C% — 5. Another example is a planar 



external field (degenerate polyhedron with two faces - 
polygon): p = 2, C2 = 2, and J m = \n± m/2\, m = 0, 1. 
This case was studied in Q for polygon=square and in || 
for polygon=line (easy plane, easy axis anisotropy). 

There is an extra symmetry in systems with a perfect 
Hamiltonian which may relate energy levels of different 
classes or inside of a class. It is the change of sign of w. 
This transformation should inverse energy levels inside 
of each class. On the other hand, all physical quantities 
depend only upon gauge invariant combinations of a type 
w k cos(JTifc), where ilk is the solid angle subtended by a 
closed loop of k > 1 walks along the geodesies. In some 
cases (polyhedra made of even-sided faces: cube, square) 
k is always even and this leads to the symmetry of the 
levels inside of each class, namely, the levels must come 
in pairs of opposite sign ±E. In other occurrences (poly- 
hedra made of odd-sided faces: octahedron, tetrahedron) 
k may be both odd and even. Then simultaneous change 
of sign of w — > — w and shift of J — > J + p/A leaves the 
invariant combinations unchanged. Therefore, each level 
E in the class of J has its counterpart —E in the class of 
J + p/4. If J and J + p/4 are equivalent, their spectrum 
is symmetric (e.g., class 3 in Table ||). 

Let us consider the magnetic response of different 
spins. We assume that magnetic field H is small, i.e., 
<7<//is|H| <C U(J) where U(J) is the variation of the 
CEF potential. The magnetic part of the Hamiltonian is 
a diagonal matrix: 

Km = J diag(-h z , h z , -h x , h x , -h y , h y ), h a = gfj, B H a , 

where g is the gyro-magnetic ratio and a — x, y, z. The 
full Hamiltonian (H + Hm) can be diagonalized for some 
symmetric directions of the field, e. g., along an easy 
direction (0,0,1). We present the ground state energies 
up to the second order in magnetic field for all classes in 
Table 

Note that class 2 has anisotropic magnetic moment, 
and classes 1, 2, and 3 have anisotropic magnetic sus- 
ceptibility in the ground state. Some classes (2, 3, and 
4) have non-zero magnetization in the ground-state. The 
low temperature behavior of these classes should be strik- 
ingly different from those without magnetization. We will 
discuss it later. 

For the 8-fold or cubic case the easy positions of CEF, 
on the unit sphere, can be chosen in the following way: 
(± arccos( l/y/3), (2fc+l)7r/4), k = 0, 1, 2,3 and enumer- 
ated as (4k + 3 ± l)/2 (see Fig. ||). The parameter space 
can be covered with p — 6 plaquettes and we should ex- 
pect Cq — 4 classes with members J m = |3n ± m/2\, 
to = 0, 1, 2, 3 all sets of levels are symmetric. The solid 
angle subtended by each plaquet is 27r/3 which gives 6 (5 
independent) equations for the phase shifts. 
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FIG. 2. Paths of the spin on the unit sphere between the 
easy positions of the field. The case of the 8-fold coordination. 

The eigenvalues of the Hamiltonian can be expressed 
in the following closed form: 



El = ±2<(tt( J 
= (3 + 2 cos x cos 



3k)), k 
2x 

y + 



-0, 1, 2, 3, 

X X : 

4 COS -COS -/(£))" 



(5) 



f(x) = (4 sin 2 I sin 2 I 
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The ground state energies of the cubic coordination 
classes and their degeneracies are presented in Ta- 
ble III. The responses of different classes to applied 
magnetic field should be quite different. To demon- 
strate this we apply magnetic field along one of the 
easy direction of the CEF, e. g., the first position 
on Fig. | (H = #(1,1,1)/V3). The Zeeman term 
of the Hamiltonian, in the chosen parameterization, is 
— hJ{3, — 3, 1, — 1, — 1, 1, 1, — 1)/3. Taylor expansions of 
the ground state energies up to the second order in mag- 
netic field are given in Table IV. There are three classes 
with non-zero magnetization in the ground state. Class 
4 has anisotropic magnetic moment, and class 3 has 
anisotropic magnetic susceptibility in the ground state. 

There exists one more type of the coordination in the 
cubic CEF — 12-fold. The 12-fold coordination is real- 
ized if the minima of the CEF are located at the cen- 
ters of the cube's edges and can be parameterized on 
the unit sphere as (kn/4, (2k + I - l)7r/4), k = 1, 2, 3, 
I — 0, 1, 2, 3. Unfortunately, the tunneling trajectories 
do not coincide with the geodesies. As a consequence, 
the number of classes, in general, is infinite. We will not 
consider the 12-fold coordination here due to its strong 
dependence on the CEF parameters. 

Now we proceed to the temperature variation of the 
magnetic response. For temperatures k B T > w large 
moments become purely classical and exhibit Curie mag- 
netic susceptibility: 

xc = (MB. 9 J) 2 /(3fc B r). 



(6) 



For temperatures k B T < w the quantum effects change 
the response drastically. The susceptibility saturates for 
the classes without magnetic moment in the ground state 
to the value: 



Xs = — 



9o 



i=l 



(7) 



where go is the degeneracy of the ground state and Xi is 
the susceptibility of the ith member of the ground-state 
multiplet. For the classes with the magnetic moment in 
the ground state Curie susceptibility persists at k B T < 
id, but its slope is different: 



Xi = —y^rni/iknT), 



(8) 



where m, is the moment of the ith member of the ground- 
state multiplet. We collect the low temperature suscep- 
tibilities of the 6- and 8-fold systems in the last column 
of Tables || and IV respectively (a common factor is 
(gJ/iB) 2 )- The exact dependencies of the susceptibility 
versus inverse temperature are given on Fig. ^ and |] for 
the 6- and 8-fold coordinations respectively. The suscep- 
tibility is isotropic in all cases due to the cubic symmetry. 
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FIG. 3. Susceptibility versus inverse temperature in the 
6-fold coordination 
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FIG. 4. Susceptibility versus inverse temperature in the 
8-fold coordination 

Now, let us deliberate experimentally observable dif- 
ferences among the classes. First of all we need to have a 
rare-earth compound with the octahedral coordination. 
The general form of the CEF potential of the octahedral 
coordination, as it was shown in ||, reads: 



U(J X , Jy, J z) — A(J X + Jy + J z 



|.7V 



3 



B{J X + Jy + J Z 



30 J X Jy J Z 



(9) 



where A and B are CEF constants. Higher order invari- 
ants are irrelevant for the most interesting situation of 
the rare-earth (R) ions with 4/-electrons. 

A general analysis of the CEF potential (^) shows that 
the plane of parameters (A-B) is divided into three sec- 
tors by the regions of stability of the 6-, 8-, and 12-fold 
coordinations: 



6-fold (r c ) 
8-fold (t%) 
12-fold (It) 



A < B/3, A < -3B/2, 
A > B/3, A > 35B/6, 
- 3B/2 < A < 355/6. 



where we enclose in parentheses the Bravais lattices of the 
cubic system which realize the respective coordinations. 
An example of the 6-fold coordination is RSb rock-salt 
structured family of compounds. 

The experimental observation of the discussed effects 
can be done on dilute alloys of the type R^R^^Sb, where 
R'=La, Lu (non-magnetic), and R=Tb (J = 6, class 5), 
Dy (J = 15/2, class 2), Ho (J = 8, class 1) or any 
other equivalent set of elements given in Table 0. Besides 
magnetization and susceptibility measurements discussed 
above one can study the EPR in these alloys. The spec- 
tral measurements can display the universal ratios of the 
resonance frequencies and check the values of the mag- 
netic momenta predicted by the theory. Our preliminary 
estimates show that Jahn- Teller effect is not essential in 
the conditions of the experiment proposed in this article 
at least at temperatures higher than 0.1K. The analysis 
of Jahn- Teller effect will be done elsewhere. 

This work was supported by the NSF grant DMR- 
9705812. We thank P.C. Canfield who indicated us RSb 
compounds' family. 
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TABLE II. The ground state energies in magnetic field 
along (0, 0, 1) direction and the low temperature magnetic sus- 
ceptibilities of the 6-fold coordination classes ((3 = 1/(&bT)) 



(10) Class Ground state energy 



Susceptibility 



1 -2w - J' z ti z /(3w) V(3w) 

2 -V2w-2Jh/3-V2J 2 h 2 /(27w) 2/3/9 

3 -2w - Jh/2 - J 2 h 2 /(Ww) (3/6 

4 -2y/2w - Jh/3 - V2J 2 h 2 /(27w) (3/9 

5 -4m - J 2 h 2 /(\2w) l/(6w) 



TABLE III. The ground state energies and their degenera- 
cies of the 8-fold coordination for w > (n = 1, 2, 3, . . .) 

Class Moment (J) Ground state First excited 

energy degeneracy state energy 



3n 

3n± 1/2 
3n± 1 
3n±3/2 



-3w 1 

— s/Ew 2 

-2w 3 

-V3w 4 



—w 





s/3w 



TABLE IV. The ground state energies in magnetic field 
along (1, 1, 1) direction and the low temperature magnetic 
susceptibilities of the 8-fold coordination classes 

Class Ground state energy Susceptibility 

-3w- J 2 ti 2 /(6w) VOM 

-V^w - Jh/3 - V6J 2 h 2 /(27w) (3/9 

-2w ~ Jh/2 - 13J 2 h 2 /(lUw) (3/6 

-V3w-2Jh/3-V3J 2 h 2 /{54w) 2(3/9 
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TABLE V. Distribution of R 3+ magnetic ions among the 
classes of the 6- and 8-fold coordinations. 



Class 


1 


2 


3 


4 


5 


6-fold 


Pm Ho 


Ce Nd Sm Gd 
Dy Er Yb 






Pr Eu 
Tb Tm 


8-fold 


Pr Eu 
Tb Tm 


Ce Gd Yb 


Pm Ho 


Nd Sm 
Dy Er 





TABLE I. The ground state energies and their degeneracies 
of the 6-fold coordination for w > (n = 1, 2, 3, . . .) 

Class Moment (J) Ground state First excited 



4 



